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SUMMARY 


An analysis is presented •vdiich indicates that the statistical 
theory of extreme values is applicable to the problems of predicting 
the frequency of encountering the larger gust loads and gust velocities 
for both specific test conditions as veil as commercial transport 
operations. The extreme -value theory provides an analytic form for the 
distributions of maximum values, of gust load and velocity. Methods of 
fitting the distribution are given along with a method of estimating 
the reliability of the predictions. 

The theory of extreme values is applied to available load data 
from commercial transport operations. The results indicate that the 
estimates of the frequency of encountering the larger loads are more 
consistent vith the data and more reliable than those obtained in 
previous analyses. 


INTRODUCTION 


In the investigation of the loads imposed on airplanes duo to 
gusty air, flight measurements of loads and gust velocities have been 
used to obtain estimates of the expected load experience under various 
operating conditions. For the problem of designing airplanes for 
ultimate static failure, the larger gust velocities are of particular 
interest, as they are likely to cause structural damage. Because the 
larger gust velocities are infrequent, measurements available from 
limited samples of data will generally not extend to the larger and 
critical values of gust velocity and load. Consequently, an important 
problem in these investigations is the development of techniques for 
estimating the probabilities of encountering these larger values. 

Previous analyses of data obtained from gust -load investigations 
have utilized the statistical approach and considered the measurements 
obtained as random samples for the conditions studied. In the analysis 
of records obtained from investigations with V-G recorders 
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(references 1, 2, and 3), the largest values of the significant varia- 
bles have been selected on some convenient and consistent basis. Fre- 
quency distributions of the largest values have been represented by 
fitted Pearson type HI probability curves. Experience vith the use of 
curves of this type has indicated, however, that the estimates of the 
probabilities of exceeding the larger and extrapolated values of the 
variables cannot be considered reliable. In some cases, there has 
appeared evidence that the estimated probabilities of exceeding the 
larger values were too low.' A further difficulty, resulting from the 
arbitrary selection of the curve type, has been the inability to derive 
satisfactory methods of measuring the reliability of extrapolated 
predictions . 

Recent developments in the statistical theory of extreme values 
(references 4 to 10) have indicated a somewhat more rational approach 
to the problem of predicting the probability of occurrence of the 
extreme values. In some cases, this new approach provides a more satis- 
factory solution to the statistical problem of the determination of the 
form of the distributions of largest values. In addition, methods are 
available for obtaining measures of the reliability of predicted values. 
The present paper summarizes some of these findings, indicates the 
method of application, and evaluates their applicability to certain 
gust -load problems. 


SYMBOLS 


m number of observation in order from smallest to largest 

N number of observations from which a mw-rlmnm is selected 

n total number of observations of wsYlminn values 

TJ e effective gust velocity, feet per second 

U-t true gust velocity, feet per second 

Ac normal-acceleration increment, g units 

V airspeed, miles per hour 

VI design level cruising speed, miles per hour 

V Q airspeed at which nra-rlmum acceleration increment occurs 

on V-G record 

x random variable 
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% 

*n-l 

% 

v(x) 

V(x) 


mth value from lovest of n values of 1 
second largest value of n values of x 
largest value of n values of x 
prpbability distribution of x 

cumulative probability distribution of x defined 
by f v(x)d x 


F(x) = 1 - W(x) 

V*(x) probability distribution of Tna-riTmim values of x 

V*(x) cumulative probability distribution of maxlmim values of x 

F*(x) = 1 - W*(x) 

T(x) recurrence period of x, average number of observations 

required to equal or exceed given values of x 

y reduced variable defined by y = o(x - u) 

y m mth value of n observations of y 

v*(y) probability distribution of reduced variable y 

V*(y) cumulative probability distribution of reduced variable y 

(a \}n) m reduced standard error of mth of n observations defined 
by equation ( 13 ) 


s 


sample estimate of a given by 



Sn = 


l.lfr 

a 


Sn-l 




k 

Sm 


u 


a 
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standard deviation of mth value of n observations 

- 

statistical parameter of distribution of extreme values 
defined by equation (5) 

statistical parameter of distribution of extreme 

Talues (r^wfa 

P 00 

universe mean value of random variable given by / x v(x)dx 


standard deviation of random variable defined by 


r 


(x - v(x)dx 


1/2 


t average number of hours per V-G record 


Subscripts: 

LLF design limit load factor 

Tnav TnflYiTniim value of random variable for specified sample 

^ tilde, used above a statistical parameter to designate the 

true universe value ' 

A bar is used over symbols to indicate the mean value of n values 

£x 

of a random variable given by — . 

BASIC STATISTICAL CONSIDERATIONS 


In the application of statistical methods to gust-load problems, 
the Tn»ylmnm values of such gust variables as normal acceleration and 
effective gust velocity obtained under given conditions have frequently 
been selected for analysis. The frequency distribution of the m a ximu m 
values obtained from successive samples has been utilized to obtain 
estimates of the probability of encountering the extreme values. The 
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problem of obtaining estimates of the probability of encountering the 
extreme values -would be considerably simplified by the a priori deter- 
mination of the underlying distribution of the TnaTiTnum values of the 
variables o The knowledge of the underlying distribution should also 
greatly increase the accuracy and reliability of the required estimates . 

Consideration of recent -work in the theory of extreme values has 
indicated that this theory would appear applicable to gust-load problems 
if the initial distribution of the variables satisfied certain require- 
ments. The theory of extreme values indicates that for certain initial 
distributions a limiting form exists for the distributions of the 
maximum values. In particular, for initial distributions of the 
exponential type, the limiting form of the distribution of maximum 
values has been found (reference 5 ) to be a simple analytic function. 
This case appears applicable to certain phases of gust statistics and 
will, therefore, be discussed in some detail. 

Consider a random variable x having a probability distribu- 
tion w(x) . The variable x is assumed to have no upper limit; that 
is, w(x) is greater than zero for all increasing values of x. The 
cumulative probability distribution from below is then given by 


W(x) = 



( 1 ) 


and indicates the probability that a measurement is less than a given 

value of x. Similarly, the cumulative probability distribution from 
above is defined by 


F(x) = / w(x)dx ( 2 ) 

Vi 

and indicates the probability that a measurement is greater than the 
given value of x. The probability that a value x is a maximum of N 
observations may be obtained from the product probabilities and is 
given by. 


w(x)dx ( 3 ) 


A simple analytical expression for W*(x) is of course not 
possible, in general, inasmuch as it depends on the form of the initial 
distribution. If N is very large, asymptotic solutions are possible 
for certain forms of the initial distribution w(x) . In particular. 
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if w(x) is of the exponential form, it has been shovn (reference 5) 
that: 

. W(x) = o-»- a(l ’ u) (4) 

vhere u is the expected largest value defined by 


V(u) = 1-1 


(5) 


and a is defined by . 


w(u) 

a = i - wfey 


( 6 ) 


The significance of the parameter a is indicated subsequently • 

The distribution of the maxi mum values of x may be obtained from 
equation (4) and is given by 


v*(x) = Av*(x) = oe- a ( x - u )- e ' a(l ‘ u) (7) 

dx 

The distribution of maximum values given by equation ( 7 ) is a two- 
parameter family of curves with parameters u and a.. The distribution 
has positive skewness with a modal value less than the mean. 

In order to simplify the foregoing results, the linear transforma- 
tion, defined by 


y = a(x - u) ( 8 ) 

\ 

is utilized to transform equation ( 7 ) to a simpler form. From equa- 
tions ( 4 ), ( 7 ), and ( 8 ), the following equations are obtained; 

V*(y) = W*(x) = e-e~y (9) 


and 


v*(y) = e”7" e ^ (10) 

The distributions of equations ( 9 ) and (10) are shown in figures 1 and 2, 
respectively. It is noted from figure 2 that the distribution of 
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equation (10) has one and only one ma xi mum at y = 0 and, therefore 
the mode of the distribution of largest values given by equation (7)* 
must be equal to u. 

. a known distribution, the rate of increase of the maximum 
value of a variable with increasing sample size is specified inflam^ 
as the average number of observations required in order to equal or 
exceed given values of the variable may be determined. If a func- 
tion T(x) is defined as the average number of maximum values required 
in order to equal or exceed given values of x, analysis has indicated 
that the value of x exceeded in T(x) observations will increase as 
a linear function of logg T(x). The results obtained in reference 6, 
page 178, indicate that for large values of x 

log e T(x) ~ a(x - u) ( 11 ) 


or 

loglO T(x) Z 0.43^3 a(x - u) (22) 

For values of x such that T(x) >10, this relationship is a suffi- 
ciently close approximation for most purposes . This result indicates 
that for large values of x the average number of samples necessary to 
equal or exceed given values of x converges toward a simple exponen- 
tial function of x. From, equation (ll), it is also apparent that the 
derivative of x with respect to the logg T(x) is a constant and is 
given by 1/a. The parameter a, thereby, specifies the rate of increase 
of the maximu m value of the variable with increasing sample size. 


Reliability of Prediction 

Equation (7) gives the probability distribution function of maxima 
values obtained from samples of a random variable with a specified type 
of initial distribution. A sample of m a xi mu m values may then be used to 
obtain estimates of parameters of the population of maximum values. It 
would be expected, however, that successive samples from the same 
population will indicate some differences in parameter estimates due to 
chance, which depend on the sample size. In practice, the parameters 
of the population are seldom known and sample estimates must be used. 
Methods have been developed for measuring the reliability of sample 
estimates by indicating a range within which, for a given probability 
level, the population value can be expected to lie. (See, for example 
reference 11.) * 

A simple and rapid method of indicating the reliability of sample 
estimates of the extreme -value distribution function called "control 
curves have been derived by Gumbel in references 8 and 9. These 
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control curves provide a measure, for a given probability level, of the 
range within which the true population value may be expected. Kimball 
in reference 10 has suggested a somewhat more precise, though consider- 
ablv more cornu lex, method of obtaining confidence limits for the distri- 
bution of extreme values. In the interest of simplicity, the present 
discussion will be restricted to the method given by uumbel. 

The method of obtaining control curves derived by Gumbel depends 
essentially on the determination of control intervals at discrete 
points along the distribution by the utilization of the properties of 
the distribution of the mth values. Within the range 0.1 < W*(x) < 0.9, 
the distribution of the mth values of a continuous variable having 
moderate skew distribution and possessing a mode is asymptotically 
normal around the most probable mth value (reference 5) with a standard 
error (cJ\/h) m given by; 


f ,-v l/[w*(*m)] [l - W*(*m)] (13) 

M,n " v*(V> 

t 

The function ( a vh. ) m is hereinafter called the "reduced standard error 
and is shown as a function of F*(*m) and y in figure 3- A hori- 
zontal interval around the true mth value given by 

• %- S m < % < % + S m • 


where ^ is defined by 


(o\/n)m 

aifn 


(15) 


gives a probability of about two-thirds that the mth value in a sample 
of n observed ma ximum values will lie in the enclosed interval. 

Similarly, the interval :% ± 2Sm has a probability of O.95 of 
enclosing the mth value from a sample of n observed maximum values. 

In order to obtain a control interval at the larger values, a 
fundamental property of the distribution of maximum values may be • 
utilized. The most probable largest value of n maximum values has 
been shown (reference 4) to have a distribution of the same form as the 
distribution of the n maximum values with the origin shifted to the 
right by a distance equal to the loge n and with the Standard deviation 
unchanged. Consideration of the probability distributions of extreme 
values indicates that an interval around the most probable largest 
value having a probability of about two-thirds is given by 
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Xn-S n <2h<sn + S n 

•where 


( 16 ) 


Sn = 


1.14 

a 


(17) 


Similarly, for a probability level of 0-95, the control interval around 
the most probable largest value is given by . ± ^»97 Since the 

control interval around the most probable largest value does not depend 
upon the number of observations n, Gumbel has suggested that a control 
intervai of equal width may be extended along the extrapolated portion 


In order to fill the gap in the control intervals for those values 
not near the median nor the extreme value, the following equation was 
obtained by Gumbel in reference 9 , page 11, for the control interval 
around the most probable penultimate value: 


( 18 ) 


(19) 


Vl - Sn.j < %_! < %-l + Sn_i 
where Sn_i is an approximation given by 

a, , - 0-754 

K-i) 

The control curves defined by equations (18) and (19) have a probability 

* or ° f robablll *r level of 0.95, the numerator of 
equation (19) becomes 1-73 instead of the value of 0 . 754 . 

Inasmuch as the foregoing method allows the deteimination of 
control intervals around the sample distribution at discrete points 
along the distribution, fairing a smooth curve through the points is 
necessary to obtain a continuous control curve. 

Closely related to the control intervals is the problem of the 
significance of observed differences between the distributions of 
samples of data. In the analysis of gust -load data, it is frequently 
required to determine whether the differences in probabilities of 
occurrence of extreme or critical values between two samples are* 
significant. The following procedure based on the control curves 
presented herein is suggested as a test of significant differences: 

(1) Obtain the control intervals with a probability of 0-95 for the 
two samples at the value of x at which the comparison is to be made. 
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(2) If each of the control intervals does not enclose the prob- 
ability value obtained from the other sample, the observed differences 
may be considered significant. 

(3) If either or both control intervals do enclose the value 
obtained from the other sample, the differences cannot be considered 
real and may be attributed to sampling fluctuations. 

Although the foregoing test for significant differences cannot be 
considered rigorous, it. would appear to be a reasonable test designed in 
accordance with the levels of significance commonly used in statistical 
tests. 


Determination of Parameters 

The preceding analysis has presented a basic distribution for 
fitting distributions of maximum values provided that the initial dis- 
tribution is of the simple exponential type. The process of fitting 
the distribution of. extreme values, to observed frequency distributions 
requires the estimation of the parameters u and a from the sample 
data. Several methods are available for the determination of the 
parameters. The accuracy with which the parameters can be estimated 
depends upon the number and accuracy of observations available. In 
practice, a minimum of about twenty-five observations has been found 
generally to be required. 

For large samples of data, n > 75, the relations obtained for the 
parameters from the method of moments have yielded satisfactory results. 
The method of moments gives the following asymptotic relations for the 
required parameters (reference 5) : 



( 20 ) 


I = 4k ^ (21) 

a it 

where c is Euler's number and equals 0.5772 and p and a are the 
mean and standard deviation, respectively, of the universe distribution 
of maxim um values. If equations (20) and (21) are assumed to be true 
-when the parameter values are replaced by their respective sample 
estimates where is estimated by x and a is estimated by s, the 
following equations are obtained: 

u = x - — (22) 

a 
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and 


For samples of data, estimates 
determined from equations (22) 
obtained for u and a are then utilized with the transformation 
equation (8) . For given values of x, equation (8) gives "equivalent 
values of y (values of y having the same probability of being 
exceeded as the given value of x) . The probabilities that a value 
of y and the associated x value will be exceeded may then be 
obtained simply from table X. 

For smaller samples, more precise methods of estimating the 
parameters are generally required in order to obtain satisfactory 
representation of the observed cumulative frequency distribution. One 
such method involves the utilization of the transformation equation (8). 
Observed values of x may be transformed to values of y by equating 
the observed points of the cumulative frequency distribution with the 
distribution of F*(y). The principle of least squares may then be 
applied to arrive at estimates of the required parameters in equation (8) 
The two normal equations that are obtained from the principle of least 
squares are 


1 \Z6s 


(23) 


of the values of u and a may be 
and (23) . The values of the narameters 


Ex 

n 


u + 


na 


(24) 


Ss: = 

n 


El Ey 2 

u — + — 


na 


(25) 


and may be solved simultaneously. Gumbel (reference 9) has suggested a 
simplification of the calculation by obtaining from equations (24) 
and (25) the following equations, in terms of the sample estimates. 


x = u + ^ 
a 


( 26 ) 


1 ( B x) 

a2 (v) 


2 

2 


(27) 


For very large samples, y approaches 0-5772, Euler's number, and (s v ) 2 
approaches ir /6j these values are the same as those obtained from the 
method of moments . 
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In the analysis of small samples, the procedure to be followed in 
determining equivalent values of y for each of the observed values 
of x requires the enumeration of the observations in order of size 
from smallest to largest. For given values of x, the ratio 

of n ~ m +' 1 gives the proportion of observations equal to or greater 
n 

than the given value. This ratio, sometimes called the recurrence 
ratio, gives a measure of the probability F*(x) that the given value 
of x will be eq ual ed or exceeded. Table I may be used simply for the 
conversion of observed recurrence ratios as observed estimates of F*(y). 
One observation is lost by this procedure because the recurrence ratio 
for the smallest observation is equal to 1 and an equivalent value of y 
cannot be determined. The loss of the smallest value is of little 
consequence because the .principal interest lies with the other end of 
the distribution. 

o 

The values of the parameters obtained by the use of equations (26) 
and (27) can then be used with the transformation equation (8) to 
obtain the probability distribution in the same manner as previously 
indicated for the case of large samples. 

In order to illustrate the application of the foregoing methods 
to gust -load data, two examples are selected and the methods of calcula- 
tion indicated for the case of both the large and small samples. 

Example 1 .- Gust -velocity measurements were available for 485 trav- 
erses of thunderstorms from the 1946 operations of the United States 
Thunderstorm Project. The frequency distribution of the maximum values 
of gust velocity per traverse is shown in table II. ' The relations 
obtained from the method of moments, equations (22) and (23), are used 
to dete rmin e the values of the parameters u and a as indicated in 
the table. By utilizing the values of the parameters, the transforma- 
tion equation given by 

u e - 12.8370 

7 4^8263 

is obtained. For given values of U e , transformed or equivalent values 
of y are obtained. Table II also indicates the equivalent values 
of y for given values of U e from 2 to 48 feet per second. The 
associated probabilities of exceeding the indicated values of U e , . 
obtained from table I, are also shown. Figure 4 illustrates the fitted 
extreme -value probability distribution along with the cumulative -data 
points . 
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Control intervals for 68-percent and 95-percent probability levels 
vere determined by means of the relations given by equations (13) 
to (19) and the results obtained are shovn in table III. The control 
vere then used to obtain the faired control curves shovn in 

figure 4. 

Example 2.- Records obtained during recent operations of a modern 
transport airplane vere selected for the second example. Tventy-six 
V-G records, each representing roughly 250 hours of flight operation 
vere available for a particular airline operator and route • The twenty- 
six m a xi m u m accelerations obtained from the records are given in 
table IV. The necessary operations in the evaluation of the param- 
eters u and a are also indicated in the table. The values of the 
parameters obtained give the folloving relation for the transformation 
equation: 


_ An - 0.8674 
7 ~ 0.2874 

The transformation equation and table I vere then used to obtain the 
probabilities of exceeding given values of An, and the results obtained 
are shown in the table. The probabilities of exceeding given values 
of An along with the cumulative data points are shown in figure 5. 

The relations for the control intervals given by equations (13) 
to (19)- were used to obtain the 68-percent and 95 -percent control 
intervals for several values of F*(An) and the results obtained are 
shown in table V. The continuous control curves obtained by fairinst 
are. shovn in figure 5. 


APPLICATION AND RESULTS ' - 


The fundamental variables generally considered in gust-load investi- 
gations include: the effective gu3t velocity U e , the normal- 

acceleration increment An, the true gust velocity U t , and airspeed V. 
In order to determine whether the distribution of extreme values is 
applicable to distributions of maximum values, consideration of the 
initial distributions of the variables is required. The available data 

on the distribution of the aforementioned variables have consequently 
been examined. 

Distributions of effective gust velocity obtained under a variety 
of operating conditions have been reported in reference J 2 . At the 
present time, however, the only extensive distributions of gust 
variables have been obtained in cloud flight from the operations of the 
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U.S. Weather Bureau Thunderstorm Project (reference 13) and the NACA 
airplane investigations of reference 14. The relative frequency dis- 
tributions of U e , An, and Ut obtained from these investigations are 
presented in figures 6, 7 , and 8, respectively. 

Extreme -value distributions were fitted to distributions of m a ximu m 
values of An available from V-G records, in order to compare the 
results obtained by the use of the extreme -value theory with those 
obtained in previous analysis of V-G data. Thirteen sets of records 
were available for analysis although six of these had only 11 
to 19 records each. The remaining seven sets had more than 25 records 
with one exception, for which only 24 records were available. One of 
these seven sets was discarded because during one flight a hurricane 
was encountered which resulted in an unusually large acceleration. This 
unusual occurrence appears to preclude a simple statistical treatment 
of the data for this set. 

Extreme-value distributions were fitted to the six sets of accelera- 
tion increments in accordance with the methods previously outlined. The 
average miles of flight required to exceed given values of acceleration 
increment were determined by. using the relation 

. O.SVr 

Mileage = — t 

■where P is the probability that a given maximum value for a record 
will exceed the given value, t is the average flight hours per record, 
and the airplane is assumed to fly at an average speed of 0.87^. 

As a basis of evaluation of the results obtained, the total miles 
of flight represented by each set of data are compared with the estimated 
number of miles (obtained from the probability curves) required to 
exceed the largest value of acceleration increment actually measured in 
each set. Figure 9(a) presents the results obtained for the six sets 
of V-G data. The straight line shown in the figure, indicates equality 
between predicted and actual mileage. For purposes of comparison, a 
similar plot is shown ’in figure 9(b) for the same sets of data but with 
the predicted mileage values obtained by means of fitting Pearson 
type III probability curves to the distributions of maximum acceleration 
increments . 

Figure 10 presents a comparison of the average life to limit load 
factor as obtained by using extreme -value and Pearson type III distribu- 
tions for the six sets of data. Inasmuch as the estimates are generally 
extrapolations, a measure of the degree of extrapolation, the ratio of 
the maximu m acceleration increment to the limit acceleration increment, 
is also given for each sample in the figure. 
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In addition to the application of the extreme -value distribution 
to distributions of maximum values of Ztti from V-G records, some 
efforts have been made to apply this distribution to the observed 
distributions of maximum values of U e and Vmnv obtained from 
V-G records. Although the results obtained appear indicative, they do 
not appear to warrant presentation in detail at this time. 


DISCUSSION 


The foregoing results have indicated that the distribution of 
extreme values given by equation ( 7 ) is the limi ting form for distribu- 
tions of maximum values where the maximum values are selected from large 
samples from an initial distribution of the exponential type. Methods 
have been presented for the fitting of extreme -value distributions. The 
methods require the estimation of only the first two moments of the dis- 
tribution and may be applied simply and rapidly . Methods have also been 
presented which allow the estimation of the reliability of the predicted 
probabilities of exceeding the extreme values. It would appear, there- 
fore, that, if applicable, the use of distribution of extreme values for 
analysis of gust loads would offer significant advantages. 

In order to determine whether the distribution of extreme values is 
applicable to gust-load variables, the available distributions have been 
examined. Distributions of Ug are available for a wide variety of 
test conditions covering turbuleiice within clear air, stratus clo uds , 
cumulus clouds, thunderstorm clouds, and areas of radar echo. The 
available data, although admittedly limited, nevertheless indicate that 
for each of the test conditions the distribution of effective gust 
velocity appears of the exponential type -with variations in average 
spacing between gusts and variations in the standard deviations for the 
different distributions. Consideration of the most extensive distribu- 
tions of U e shown in figure 6 indicates that the exponential dis- 
tribution shown by the fitted straight line are good representations 
of the data. Therefore, the distribution of extreme values would appear 
to apply bo the effective gust velocity for given test conditions. 

The available distributions of An have been obtained, largely, 
under test flight conditions in which the airspeeds have been kept 
constant or at least restricted. The linear relation between U e 
and An in the sharp -edge gust relation would, consequently, be 
expected to yield the same form for the distribution of U e and An. 

The distribution of An shown in figure 7 substantiate the expectation 
that the form of the An distribution is the same as the U e distribu- ' 
tions of figure 6 . It may be concluded that for test flights at 
constant speed at least the distribution of An would remain in the 
exponential form and the distribution of extreme values would be 
applicable. 
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The available evidence as indicated by figure 8 appears to indicate 
that the distribution of U-fc also follows .a simple exponential distribu- 
tion. The distribution of extreme values may, therefore, also be applied 
to distributions of maximum values of U-fc. 

The question of the applicability of the extreme -value distribution 
to the maximu m values of U e obtained from V-G records from commercial 
operations arises since these records generally cover a wide variety of 
operating and weather conditions. Although the available data, in some 
cases, appear to indicate that the distributions of U e are of the 
exponential form, they do not extend to sufficiently high, values to be 
conclusive. Available information on commercial -transport gust, exper- 
ience and on the frequency distributions of U e for various weather 
conditions indicates that the larger values encountered in commerical 
operations result largely from flight through thunderstorms and con- 
vective clouds. On the basis, then, that the distribution of gusts in 
convective clouds is exponential, it seems reasonable to expect that 
the distribution of U e in commercial operations should have an exponen- 
tial distribution at least for the larger values of gust velocity. 
Application of the distribution of extreme values would seem reasonable, 
therefore, to the distribution of maximum values of U e obtained 
from V-G records. 

In connection with the possible application of the distribution of 
extreme values to maximum values of acceleration increment obtained 
from V-G records of commercial transport operations, the question of the 
airspeed operating practice may be of some importance. Apparently, a 
systematic variation of airspeed as a function of gustiness could con- 
ceivably cause an appreciable departure from the exponential form in the 
distribution of even' at the larger values. Little evidence is 

available, however, to indicate that such a systematic variation exists 
in practice . On the assumption, then, that operating speeds in areas of 
moderate to severe turbulence are largely confined to a narrow speed 
range for a given airplane or are independent of the intensity of the s 
turbulence, the application of the distribution to the maxi mi rm values of 
acceleration increment is reasonable. 

Consideration has also been given to the applicability of the 
extreme -value distribution to maximum values of airspeed obtained from 
V-G records. Little information is available, however, concerning the 
speed-time distribution of airspeed. Several attempts were made, con- 
sequently, to fit observed distributions by using the distribution of 
extreme values. The results have indicated that the observed distribu- 
tions of Vtopy are frequently skewed negatively as compared with the 
positive skew of the distribution of extreme values. As a result, the 
agreement between, the observed data and the distribution of extreme 
values was poor. The limited information available on the speed-time 
distribution of airspeed has further indicated that the distribution of ■ 
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airspeeds follows no simple exponential form and depends largely on 
operational practice. The distribution function of equation (7) would, 
therefore, appear not to he applicable to distributions of m»Tinni m 
values of airspeed. 

Another variable frequently studied in gust -load investigations is 
the airspeed at which the maximum acceleration on a V-G record is 
encountered. By definition, this variable does not represent maximum 
values of airspeed and, as a result, no initial distribution that is 
directly related to this variable is apparent. Available distributions 
of V 0 have been e x a min ed, however, and the results indicated that the 
distribution of extreme values of equation (7) is not applicable to dis- 
tributions of these data. 

It has been indicated that the distribution of extreme values of 
equation (7) is a rational distribution form for purposes of fitting 
distributions of m a ximum , values of 4n, U©, and TJ-fc when the maximums 
are selected from a large number of observations. The need remains for 
indicating that the application of this distribution yields accurate 
predictions. The verification of statistical techniques of this type 
requires huge masses of data. Sufficient data are not available to 
enable a complete answer to this question. Available data may be used, 
however, to compare the estimates obtained with those obtained by the 
methods commonly used in the analysis of V-G records. 

The comparison shown in figure 9(a) between the estimated miles to 
exceed the highest observed value of An based on extreme -value methods 
and the actual mile3 flown indicates that the agreement is extremely good 
In only one case, is any appreciable difference noted and in this case, 
the actual mileage is about twice the estimated mileage. The same com- 
parison between the estimated miles to equal or exceed the maximu m 
values obtained by Pearson type III curves and the actual miles flown is 
shown in figure 9(b) • The figure indicates appreciably more scatter 
with the predicted miles always greater than the actual miles with the 
ratio of the predicted to actual miles in two cases about 3 to 1. This 
tendency toward overestimation of miles to exceed given acceleration- 
increment values is unconservative and may be misleading. On this basis 
the extreme -value distribution apparently yields more reliable estimates* 
of the frequency of encountering the larger values. 

Since the miles required to exceed limit load factor is of 
particular concern, the comparison shown in figure 10 of the estimates 
to exceed limit load by th9 extreme -value and Pearson type I El curves 
becomes of interest . Consideration of figure 10 indicates that the 
Pearson curves in every case yield larger mileage estimates. The 
dift erences between the estimates at limit load factor vary appreciably 
as indicated by the ratio of the mileages which vary from about 1.3 
to 60. Consideration of the ratios of the maximum values of An 
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observed to limit load factor for each set of data Indicates that the 
magnitude of the differences between the results of the two methods 
increases with the degree of extrapolation required to reach limit load 
factor. 

In summary, then, the foregoing comparisons have indicated that as 
compared with previous methods of analysis, extreme -value methods yield 
results more consistent with the data at the largest values . The 
differences between the estimates are considerably amplified when com- 
parisons are made at extrapolated portions of the distributions. In 
view of this evidence and the rational foundation for the extreme-value 
distribution presented previously, the use of this distribution would 
generally be expected to provide more reliable estimates. 

In the practical application of extreme -value methods to the 
analysis of V-G records, it should be noted that the extreme-value 
distribution of the maximum values of £n may be used, in conjunction 
with the methods of reference 1, to predict the velocity-acceleration 
or V-N envelope for a given airplane. It has also been suggested that 
the determination of the V-N envelope should be made by analysis of the 
accelerations by speed bracket. In this connection, the use of extreme- 
value methods would also seem proper provided that for each record a 
large number of accelerations are experienced within each of the speed 
brackets chosen. 

Another problem encountered, in practice, is that estimates of 
the frequency distribution are frequently required for very limited 
samples, sometimes as few as ten records. In order to increase the 
samcle size, the practice has been to use both the ma x i m um positive and 
maximum negative acceleration increments from each record; the sample 
size is thereby doubled. Although this practice departs somewhat 
from the extreme-value theory, an analysis of available data indicated 
that the use of two maximums per record for the larger samples had only 
a minor effect on the estimates obtained. The use of two m a x im u m s per 
record would appear, therefore, to offer a useful basis for increasing 
sample size, particularly for very limited samples. 


CONCLUDING REMARKS 


An analysis of gust -load data by the use of the statistical theory 
of extreme values has indicated the following: 

1. The theory of extreme values gives a rational form for the dis- 
tribution of maximum values that appears applicable to distributions of 
TmnYirmim values of the effective gust velocity, true gust velocity, and 
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normal-acceleration increment obtained from test flights or commercial 
transport operations. 

2. Simple methods are available for fitting the distribution of 
extreme values to samples of data and for obtaining control curves that 
provide a measure of the reliability of the estimates of encountering 
the larger values. A simple and reasonable, but nonrigorous, method of 
determining -whether differences betveen two estimates may be considered 
significant is also presented. 

3 • The application of the distribution of itiw-vimuTn values to 
available V-G data yields estimates of the frequency of encountering 
the larger acceleration increments that are consistent -with the avail- 
able data and appear to be more reliable than the estimates obtained in 
previous analyses • 

4. Notwithstanding the fact that the methods described in the pre- 
sent paper appear to be applicable, the possibility must always be kept 
in mind that because of limitations of the method, available records, 
or both, a value may occur that far exceeds the limi ts of the 
theoretical predictions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., July 22, 1949 
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TABLE I 

CUMULATIVE PROBABILITY DISTRIBUTION F*(y) = 1 - e" e ” y 


7 

F*(y) 

7 

F*(y) 

7 

F*(y) 

7 

F*(y) 

-2.50 

0.999995 

0.15 

0.5771 

2.80 

0.05900 

5.45 

0.004287 

-2.45 

.999991 

.20 

•5590 

2.85 

.05620 

5.50 

.004079 

- 2.40 

.999984 

.25 

.5410 

2.90 

.05353 

5.55 

.003880 

-2.35 

.999972 

.30 

.5233 

2.95 

.05100 

5.60 

.003691 

-2.30 

•999953 

.35 

.5057 

3.00 

.04856 

5.65 

.003512 

-2.25 

•999924 

.40 

.4884 

3.05 

.04625 

5.70 

.003341 

-2.20 

.99988 

•45 

.4714 

3.10 

.04405 

5.75 

.003178 

-2.15 

•99981 

•50 

.4547 

3.15 

.04194 

5.80 

.003024 

-2.10 

.99972 

•55 

.4384 

3.20 

•03994 

5.85 

.002876 

-2.05 

•99958 

.60 

.4224 

3.25 

.03803 

5.90 

.002735 

-2.00 

.99938 

.65 

.4067 

3.30 

.03621 

5-95 

.002603 

-1.95 

.9991 

•70 

•3914 

3.35 

.03447 

6.00 . 

.002476 

-1.90 

.9988 

•75 

.3765 

3.40 

.03282 

6.05 

.002355 

-1.85 

.9983 

.80 

• 3619 ' 

3-45 

•03125 

6.10 

.002241 

-1.80 

.9976 

.85 

.3478 

3.50 

•02975 

6.15 

.002132 

- 1-75 

.9968 

.90 

• 3341 

3.55 

.02831 

6.20 

.002027 

-1.70 

.9958 

•95 

.3207 

3.60 

.02695 

6 .25 

.001928 

-I.65 

.9945 

1.00 

.3078 

3.65 

.02565 

6.30 

.001834 

-1.60 

.9929 

1.05 

.2952 

3.70 

.02442 

6.35 

.001746 

- 1.55 

.9910 

1.30 

.2832 

3 - 75 . 

.02325 

6.40 

.001661 

- 1.50 

• 988 r 

1.15 

.2714 

3.80 

.02212 

6.45 

.001579 

-1.45 

•9859 

1.20 

.2601 

. 3-85 

.02105 

6.50 

.001502 

- 1.40 

.9827 

1.25 

.2491 

3.90 

.02004 

6-55 

.001429 

- 1-35 

.9789 

1.30 

.2385 

3-95 

.01907 

6 .60 

.001359 

- 1.30 

.9745 

1-35 

.2283 

4.00 

.01815 

6.65 

.001293 

-1.25 

.9695 

1.40 

.2186 

4.05 

.01727 

6.70 

.001230 

-1.20 

•9638 

1.45 

.2091 

4.10 

.01643 

6 . 75 - 

.001170 

- 1.15 

•9575 

1.50 

.2000 

4.15 

.01564 

6.80 

.001113 

-1.10 

•9504 

1-55 

.1912 

4.20 

.01489 

6.85 

.001059 

- 1.05 

.9426 

1.60 

.1828 

4.25 

.01416 

6.90 

.001008 

-1.00 

•9340 

1.65' 

•1747 

4.30 

.01348 

6.95 

.000959 

-.95 

•9247 

1.70 

.1670 

4-35 

.01283 

7.00 

.000912 

-.90 

•9145 

1-75 

•1595 

4.40 

.01220 

7-05 

.000867 

-.85 

.9036 

1.80 

.1524 

4.45 

.01161 

7-10 

.000825 

-.80 

.8920 

1.85 

.1455 

4.50 

.01105 

7.15 

•.000785 

-.75 

.8796 

1.90 

•1390 

4.55 

.01051 

7-20 

.000747 

-.70 

.8665 

1-95 

.1326 

4.60 

.01000 

7.25 

.000710 

-.65 

.8527 

2.00 

.1266 

4.65 

.009516 

7-30 

.000676 

-.60 

■ .8383 

2.05 

.1208 

4.70 

.009054 

■ 7-35 

.000643 

-.55 

.8233 

2.10 

• 1153 

4.75 

.008614 

7-40 

.000611 

-.50 

.8077 

2.15 

.1100 

4.80 

.008196 

7.45 

.000581 

-.45 

.7916 

2.20 

.1049 

4.85 

.007798 

7.50 

.000553 

-.40 

.7750 

2.25 

.1000 

4.90 

.007420 

7-55 

.000526 

-•35 

.7581 

2.30 

.09543 

4.95 

.007058 

7.60 

.000500 

-.30 

.7407 

2-35 

.09096 

5.00 

.OO6716 

7.65 

.000476 

-.25 

•7231 

2.40 

.08673 

5.05 

.OO6389 

7.70 

.000453 

-.20 

.7052 

2.45 

.08266 

5.10 

.006078 

7.75 

.000431 

-.15 

.6871 

2.50 

.07880 

5-15 

.005782 

7.80 

.000410 

-.10 

.6689 

2-55 

.07511 

5.20 

.005502 

7.85 

.000390 

-.05 

.6505 

2.60 

.07158 

5.25 

.005235 

7.90 

.000371 

0 

.6321 

2.65 

.06822 

5.30 

.004980 

7*95 

•000353 

.05 

.6137 

2.70 

. .06500 

5.35 

.004737 

8.00 

.000336 

.10 

•5954 

2.75 

.06193 

5.40 

.004507 
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TABLE H 

SUMMARY OF CAICULATIONS, 1946 THUUBERSTORM PROJECT DATA 


[Example I] 


(fpe) 

Midpoint 

Frequency, 

f 

CunrulatiVe 

frequency, 

F 

Relative 

cumulative 

frequency, 

F/N 

V 

U e 2 f 

2 to 4 

3 

4 

485 

1.000 

1 2 

36 

4 to 6 

5 

11 

481 

.9918 

55 

275 

6 to 8 

7 • 

27 

470 

.9691 

189 

1,323 

8 to 10 

9 

48 

443 

.9134 

432 

3,888 

10 to 12 

11 

62 

395 

.8144 

682 

7,502 

12 to 14 

13 

58 

333 

.6866 

754 

9,802 

14 to 16 

15 

55 

275 

.5670 

825 

32,375 

16 to 18 

17 

60 

220 

.4536 

1020 • 

17,340 

18 to 20 

19 

6l 

160 

.3299 

1159 

22,021 

20 to 22 

21 

36 

99 

.2041 

756 

15,876 

22 to 24 

23 

17 

63 

.1299 

391 

8,993 

24 to 26 

25 

18 

46 

.09485 

450 

11,250 

26 to 28 

27 

8 

28 

•05773 

216 

5,832 

28 to 30 

29 

7 

20 

.04124 

203 

5,887 

30 to 32 

31 

6 

13 

.02680 

186 

5,766 

32 to 34 

33 

3 

7 

.01443 

99 

3,267 

34 to 36 

35 

1 

4 

.008247 

35 

1,225 

36 to 38 

37 

2 

3 

.006186 

74 

2,738 

38 to 40 

39 

1 

1 

.002062 

39 

1,521 

z 

— 

485 



7577 

136,917 

*e 

(fps) 

y 

F*(y) 

Derivation of transformation equation: 

2 

-2.2454 

0 .9999 

u e = 15 

.6227 



4 

-1.8310 

.9980 




6 

-1.4166 

.9838 





8 

-1.0022 

.9344 





30 

-.5878 

.8346 





12 

-.1734 

.6956 

1 - 4.®>6"? 



14 

.2410 

.5442 

a 




16 

.6554 

.4051 





18 

I.O698 

.2904 





20 

1.4842 

.2029 





22 

I.8986 

.1392 


.1899 



24 

2.3130 

.09427 

3U e = 6 



26 

2 .7273 

.06332 




28 

3-1417 

.04229 




* 

30 

3-5561 

.02814 





32 

3-9705 

.OI869 

u = 12 .8370 

^ MACAWS '' 

34 

4-3849 

.01239 



36 

4-7993 

.008202 





38 

5-2137 

.005429 





40 

5-6281 

.003590 





42 

6.0425 

.002373 


12.8370 



44 

6.4569 

.001568 

h 



46 

6.8713 

.001037 



% 

48 

7.2857 

.000687 
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TABLE IV 


SUMMARY OF CAICULATIONS, V-G DATA 


[Example 2] 


An 

(g units) 

Serial 

number, 

m 

n - m + 1 

7 

y 2 

Ar|2 

n 

0.50 

1 

1.0000 



0.2500 

.61 

2 

.9615 

-1.182 

1-3971 

.3721 

.65 

3 

•9231 

-.942 

.8874 

.4225 

.66 

4 

.8846 

. -.770 

.5929 

.4356 

.70 

5 

.8462 

-.627 

•3931 

.4900 

-73 

6 

.8077 

-.500 

.2500 

.5329 

.80 

7 

.7692 

-.383 

.1467 

.6400 

.80 

8 

•7308 

-.272 

.0740 

.6400 

.82 

9 ■ 

•6923 

-.164 

.0269 

.6724 

•91 

10 

•6538 

-.059 

.0035 

.8281 

•91 

11 

.6154 

.045 

.0020 

.8281 

.98 

12 

.5769 

.151 

.0228 

.9604 

1.00 

13 

.5385 

.257 

• .0660 

1.0000 

1.00 

14 

.5000 

.366 

.1340 

1.0000 

1.00 

15 

-4615 

.480 

•2304 

1.0000 

1.05 

16 

.4231 

•598 

•3576 

1.1025 

1.05 

17 

•3846 

.723 

.5227 

1.1025 

1.10 

18 

•3462 

.856 

.7327 

1.2100 

1.20 

19 

•3077 

1.000 

1.0000 

1.4400 

1.25 

20 

.2692 

1.160 

1.3456 

1.5625 

1.30 

21 , 

.2308 

1.338 

1.7902 

1.6900 

1.30 

22 

.1923 

1.544 

2.3839 

1.6900 - 

1*35 

23 

•1538 

1.790 

3 .2041 

1.8225 

1.1*0 

24 

.1154 

2.099 

4 .4058 

1.9600 

1.65 

25 

.07692 

2.525 

6.3756 

2.7225 


1 .oO 

26 

.03846 

3.239 

10.4911 

3.2400 

X 




13 .272 

36.8361 

29.6146 1 







J 


An 

7 

F*y 

0.50 

-1.2784 

O.9723 

.60 

-.9304 

.9207 

.70 

-.5825 

.8331 

.80 

-.2345 

.7176 

.90 

.1134 

.5905 

1.00 

* .4614 

.4676 

1.10 

.8093 

•3593 

1.20 

1.1573 

.2698 

1.30 

1-5052 

.1991 

1.40 

1.8532 

.1451 

1.50 

2 .2011 

.1048 • 

1.60 

2 .5491 

.07518 

1.70 

2.8970 

•05369 

1.80 

3 .2450 

.03822 

1.90 

3.5929 

.02714 

2.00 

3-9408 

.OI925 

2.10 

4.2868 

.01363 

2.20 | 

4.6367 

.00964 


Derivation of transformation equation: 
Ai = 1.0200 

B An = 0.3202 

-i = 0.287t 

a 

y = 0.5309 

By = 1.1141 X[naca^" 

u = 0.8674 


y = 


An - 0.8674 
0.2874 















30 


NACA TN 1926 



Figure 4.— Probability that the maximum yalue of U e per traverse will 
exceed indicated value; 1946 Thunderstorm Project data; example 1. 
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Figure 5* — Probability that the maximum value of acceleration increment 
per record will exceed indicated value j V—G data] example 2. 
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(a) 1946 Thunderstorm Project data. Total number of observations, 27,486. 
Figure 6.— Distribution of effective gust velocity U e . 
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Figure 6.— Continued. 
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Figure 6.— Concluded. 
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(a) 1946 Thunderstorm Project data. Total number of observations, 21,385- 
Figure 7*— Distribution of normal--acceleration increment An. 
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(b) 19^7 Thunderstorm Project data. Total number of observations, 17*367* 


Figure 7*— Continued. 
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Art } g un/ts 

(c) 1941 and 1942 data. Total number of observations, 4645. 

Figure 7.- Concluded. 
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(a) 1946 Thunderstorm Project data. Total number of observations, 6715. 
Figure 8.- Distribution of true gust velocity U t - 
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Ui^fps 

(b) 19V7 Thunderstorm Project data. Total number of observations, l8l6. 

Figure 8.— Continued. 











NACA TN 1926 


4i 



> (a) Extreme values. 


Figure 9*— Miles flown as 
normal— ac c e ler at i on 


a function of miles to exceed, largest observed 
increment for six samples of Y-G records. 


* 











Figure 10.— Comparison of. two methods of predicting life to limit load factor from V-G data hy use of 

one maximum. 




